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Abstract 

Studies are made of the turbulent separation bubble in 
a two-dimensional semi-infinite blunt plate aligned to a 
uniform free stream when the oncoming free stream 
contains a pulsating component The discrete-vortex 
method is applied to simulate this flow situations because 
this approach is effective to represent the unsteady 
motions of turbulent shear layer and the effect of viscosity 
near the solid surface. The numerical simulation has fairly 
reasonable predictions with the experimental results which 
have already performed. A particular frequency given a 
minimum reattachment which is related to the drag 
reduction and the most effective frequency is dependent 
on the most amplified shedding frequency. The turbulent 
flow structure is scrutinized, which includes the time-mean 
and fluctuations of the velocity and the surface pressure, 
together with correlations between the fluctuating 
components. A comparison between the pulsating flow 
and the non-pulsating flow at the particular frequency of 
the minimum reattachment length of the separation bubble 
suggests that the large-scale vortical structure is associated 
with the shedding frequency and the flow instabilities. 


Introduction 

An improved understanding of pulsating flow 
characteristics on the turbulent separation bubble is useful 
in the design of aerodynamic high-speed vehicles and of 
pulsating turbomachineries. In this study, the separation 
bubble is generated by flow separation from a sharp 
leading-edge of a blunt circular cylinder whose axis is 
aligned parellel to the approaching main flow. Relevant 


experimental studies have been earned out by Kiya ct al. ! 
with a view toward the control of a turbulent leading-edge 
separation bubble. 

Control of the separation bubbles by sinusoidal 
perturbations has been reported by many researchers 2 '*. 
Roos and Kegelman 2 obtained the reattachment length in a 
backward-facing step flow as a function of the level and 
frequency of the perturbation. Sigurdson and Roshko 
analyzed the structure and control of a turbulent 
reattaching flow; the reduction of the pressure drag, 
bubble height and reattachment length were found to 
depend critically on the forcing frequency. 

In order to numerically simulate pulsating flow on a 
turbulent leading-edge separation bubble, the discrete- 
vortex model is applied. This approach has been 
demonstrated to be effective in representing the unsteady 
motions of turbulent separation bubble 4 . The finite- 
difference simulation of the averaged Navier-Stokes 
equation with turbulence models has the advantage in 
computational accuracy and applicability. However, it is 
currently limited to the ranges of low Reynolds numbers 
for separated unsteady flows. 

The discrete- vortex model is known to be a powerful 
tool for simulating unsteady separated flows of high 
Reynolds numbers 0 . The separating flow at the leading 
edge is represented by a combination of an inviscid 
potential flow and discrete vorticies. The effect of 
viscosity near the solid surface is incorporated in the 
model. The reduction in the circulation of elemental 
vortices is also introduced as a function of their ages in 
order to represent the three-dimensional deformation of 
vortex filaments. Details regarding the numerical 
procedures can be found in Ref. 5. The ability of the 
discrete vortices to adequately represent pulsating 



(a) Pulsating flow 



Acoustic speaker 


(b) Experimental apparatus with acoustic speaker 
Fig. 1 Experimental apparatus and numerical analogy 

continuous vortex sheetes in the separating bubbles has 
been tested in Ref. 6-7. 

In this paper, we consider the specific example of the 
flow past a two-dimensional flat plate with finite thickness 
and a blunt leading-edge, which is aligned parallel to a 
uniform approaching stre am. The pulsation was provided 
by a sinusoidally oscillating jet issued from a thin slit 
along the separation edge. Therefore, in the numerical 
simulation, the oncoming free stream was assumed to 
contain a sinusoidally-varying pulsating component The 
entire flow field in the separation bubble is affected 
mainly by the dynamics near the leading-edge where 
separation occurs (Fig. 1). 

The purpose of this study is to examine the details of 
the turbulent leading-edge separation bubble by discrete- 
vortex model. The two key parameters characterizing the 
free stream are the amplitude of pulsation A, and the 
frequency parameter St (=fH/Uj). The effects of the 
separation bubble on the forebody drag, the reattachment 
length and the lock-on effect between the pulsating 
frequency and the shedding frequency are investigated. 
The turbulent flow structure is also scritinized, which 
includes the time-mean and fluctuations of the velocity 


and the surface pressure, together with the correlations 
between the fluctuating components. The effect of 
pulsation on the minimum reattachment length of the 
separation bubble is examined in detail. This will show 
that the large scale vortical structure is closely linked with 
the issues of the shedding frequency and the flow 
instabilities. 

Discrete-vortex model 

The leading-edge separation bubble of a blunt two- 
dimensional body is considered. This flow geometry is 
basically the same as the flow configuration of Kiya . 
Specifics regarding the utilization of discrete-vortex 
method can be found in their studies. The separation 
bubble is generated by flow separation from a sharp 
leading-edge of a blunt two-dimensional body. The 
dimension of the flat plate is of finite thickness(H) and 
semi-infinite length. The interactions between the two 
separation bubbles at the comers arc assumed to be 
minimal. Thus, the symmetry condition is applied in this 
problem. The geometry of the body is given in Fig.2. 

The Schwartz-Christoffel transformation is used to 
project the exterior region of the body ( the physical 
plane, z-plane) into an upper half plane ( the transformed 
X-plane). The transformation is given by 

z = ■§-[(2. 2 — 1) WI — cosh 'A.]+ iH (1) 

7l 

The upstream free stream velocity U ; contains a 
pulsating component, therefore, 

1 + Asin2nft) (2) 

where f is the pulsating frequency and A is the amplitude 
of pulsation. If the flow has no pulsation, i.e., A=0, \J* 
reverts directly to the constant velocity at upstream 
infinity Uj. The complex potential w r induced by the 
discrete vortices is given by 

NilC. X-X 

w,= S^" Iog X- % (3) 

where N is the total number of vortices in the flow field, 
and j is the position of the jth vortex in the transformed 
plane with its complex conjugate and Kj denotes its 
circulation. The complex potential w for the entire flow 
field is the sum of w ; and w, by superposition, i.e., 
w=w r +w t , where w t is the complex potential for the 




(a) Physical z plane 
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where Kn is the initial strength and the position of the 
nascent vortex is assumed to be i(H+0.5e), and At, is the. 
time interval between the introduction of the nascent 
vortices. The position of the nascent vortex is assumed to 
be i(H+e), where e is an approximate initial value. The 
strength and location of the nascent vortices can be 
adjusted by an appropriate iteration scheme satisfying the 
Kutta condition 5 . 


(b) Transformed X plane 
Fig. 2 Geometry of flow 


irrotational flow around the body. 

The velocity field in the physical plane is given by 


. dw dw dX 
u “ lv= dF = dA dF 


(4) 


where u and v are the velocity components in the x- and 
y-directions, respectively. The velocity at a vortex point, 
the kth vortex, has to be obtained by differentiating with 
Taylor expansion and the complex invariant. 
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where dw^dz is deduced from taking the appropriate limit 
as z— > z^. 

The convection of vortices was advanced by a second 
order scheme with a small time interval, the nascent 
vortices have to be calculated at a much smaller time 
interval. Since the velocity at the leading-edge should be 
finite, the complex potential should satisfy the Kutta 
condition at the edge 

(dw/dX)x.i=0. 

The vortices are shed from the leading-edge into the 
separation bubble as a result of the separation. The rate of 
vorticity shedding was determined through the 
relationship, which was inferred from the experimental 
results* 


The reduction of the circulation of every vortex is 
modelled by 5 . 



where K(t) is the circulation at time t, a is an adjustable 
constant, and Re denotes the Reynolds number U,(H/v) 
and v the kinematic viscosity of the fluid. The destruction 
and coalescence of vortices are assumed to be proportional 
to the viscous core radius of each vortex. Physically, the 
removal of potential vortex may be regarded as 
corresponding to the destruction of vorticity in the shear 
layer by interaction with the boundary layers along the 
surface. The decay law (8) was deduced from the exact 
solution of Navier-Stokes equations for a single 
rectilinear viscous vortex if r is replaced by the radial 
distance from the center of the vortex. After caning out a 
number of preliminary calculation, the value of the 
product a 2 Re = 60 was employed, and this was found to 
achieve satisfactory agreement with experimental result . 
It may be noted that an optimum value of a Re depends 
on the particular type of flow considered. 

The pressure coefficient Cp can be calculated from 
the Bernoulli equation. 
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where is the pressure of the free stream, p the density 
of fluid and the velocity potential. 

In the course of computations, some vortices approach 



very closely the wall of the blunt body. Consequently, 
due to the presence of image vortices, these vortices 
would have unreasonably large velocities. In order to 
rectify this computational problem, vortices that 
approched the wall nearer than the depth of 0.02H 
were removed from the flow field. Since the transport of 
momentum and vorticity were negligible in the far field 
downstream, vortices and their images that were located 
further than the region 25H were also removed from the 
computation. 

Since a large number of vortices exist with random 
locations in the flow fleld, it is probable that some vortices 
attain small separation and, therefore, produce velocity 
jumps at each other’s positions due to the absence of 
viscosity. In order to alleviate this difficulty, the concept 
of the cut-off vortex, which was originally suggested by 
Chorin 9 , is also employed, i.e., 


K(r/o) 

2tc 


(r£o) 


( 10 ) 


where is the stream function of the vortex and c the 
cut-off radius. The numerical value c = 0.05H was 
adopted in this simulation . The concept of the cut-off 
vortices was justified in that the vortex blob in discrete- 
vortex simulation is basically different from that of 
potential flow theory, say, the point vortex. 

It is noted again that the viscosity of fluid is 
instrumental in enforcing the no- slip flow condition. 
Since the discrete-vortex simulation is started from the 
in vise id flow, an appropriate procedure should be devised 
to include the viscous effect for the turbulent separation 
bubble and the shear layer. The viscosity of fluid gives 
rise to the displacement thickness in the boundary layer 
and, thereby, transfers momentum to the direction of 
transverse velocity. For this purpose, an artificial 
transverse velocity V d (=O.0125Ui) was added uniformly 
to the edge of the shear layer. 

It seems that the vortex distributions and wave forms 
of flow field evolve to be statistically stationary at times 
tyt/H > 80. In this sense, the mean values and 
fluctuating components were calculated between the 
interval of 80 < U t t/H < 280. The non-dimensional 
time step t for the movement of the vortices was 0.16H / 
Ui and the time interval At* between the introduction of 
the nascent vortices was 0.32H / Uj ( = 2At). 



Fig. 3 Comparison of X/X* between experiment and 
computation 



Fig. 4 Reattachment length^/X*) with Strouhal 
number 



Fig. 5 Reverse flow intennittency distribution 






Results and discussion 


It is important to verify the aforementioned analogy 
between the sinusoidal oscillating jet at the separation 
edge and the pulsating ffeestream flows by the present 
discrete-vortex method. Toward this end, the reattachment 
length normalized in the form X/X ro is plotted in Fig.3 
against the non-dimensional frequency, i.e., the Strouhal 
number St for three different levels of the perturbation. It 
is noteworthy that the comparision gives a fairly 
consistent prediction, and this supports the assertion that 
this analogy seems to be reasonable. The most interesting 
feature of Fig.3 is that the reattachment length attains a 
minimum at about St=1.4 for both cases, while the 
amplitude affects only on the total size of reattachment. 
This implies that the separation bubble is affected mainly 
by the frequency of perturbation. 

In the present study, the turbulent structures are 
scrutinized with the comparisons for the non-pulsating 
flow, natural flow(A=0) and for the perturbed flow at 
A=0.3 and St=0.1. It should be noted that the condition of 
flow pulsation was selected such that a maximum 
reduction of reattachment length for this two-dimensional 
flow configuration can be realized not for the circular 
blunt body (Fig. 4). For this two cases, the reverse-flow 
intermittency is displayed in Fig. 5. The reattachment 
position is defined as the point where the intermittency, I, 
has the value of 0.5. In this figure, it is noticed that the 
reattachment position(X/H=6.5) of perturbed flow is much 
reduced in comparision to that of non-pulsating 
flow(X/H=9.1). This agrees with the experimental result 5 . 

Fig. 6 shows the distribution of surface velocity under 
the same condition. In a perturbed flow, the position of 
zero surface velocity is located slightly upstream of the 
reattachment position while the same position is 
maintain ed in natural flow. This means that the reverse 
flow intensity of perturbed flow is much stronger than that 
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Fig* 7 (a) Time-mean pressure coefficient distribution 
along the surface 



Fig. 7 (b) R.m.s. pressure distribution along the surface 

of natural flow, and the feed back of fluid toward 
upstream of perturbed flow is active. 

The distributions of pressure and its fluctuating 
pressure along the plate surface are shown in Fig. 7. This 
points to the relation with the flow reattachment for both 
perturbed and natural flows. It is known that the flow 
reattachment occurs at near the mean pressure recovery 
region. It is found that the mean pressure coefficient oi 
perturbed flow at separation line is smaller than that of 
natural flow (Fig. 7(a)). This implies that the pressure drag 
is reduced when the flow is perturbed. In the pressure 
fluctuations of Fig.7(b), this feature of reattachment 
position can be elucidated more clearly. The r.m.s. of 
pressure fluctuation is relatively small near the separation 
region and then it increases along with the downstream 
and reaches maximum just in front of reattachment 
position. This region is also the same position of the 
maximum mean pressure recovery rate. In turbulent shear 
layer, it is believed that the flow has its own natural 


Fig. 6 Surface velocity distribution 





instability and unsteadiness and the flow disturbance can 
be amplified most effectively at the appropriate pulsation. 

This is called the most effective frequency, St^. This 
effect of amplification is maximized at the flow 
reattachment region, so the motion of large-scale vortices 
in this region are the most important elements that can 
decide the characteristics of the entire flow field. 

In passing, it can be found in Fig. 7(a) that the 
pressure coefficient at separation point(Cp,) for the 
pulsating flow is reduced considerably about 30% than 
that of non-pulsating flow. Reduction of Cp, reproduces 
the reduction of pressure drag of the body. Tbe maximum 
reduction of Cp, is also observed at St^. This agrees well 
with the experimental result of Sigurdson and Roshko . 
Koenig 10 studied the relation between Cp, and the drag 
coefficient (Cp.) of fluid bodies, and it gives 

Cp.~0.840.2Cp,. 00 

From this relation, it can be noticed that the drag of the 
body (Cp.) is decreased approximately 6.25% when 
perturbed, i.e., C^,=0.68 for non-pulsating flow while 
CarKJ.64 for pulsating flow. 

The main reason of drag reduction can be explained 
in two ways. First, the pulsating perturbation invokes 
vortex coalescence in the separation bubble and then it 
causes the enhancement of flow spreading rate. Since the 
downstream flow rate is increased due to these 
phenomena, it gives the drag reduction. It means that the 
enhancement of vortex merging and flow spreading rate is 
strongly dependent upon the pulsating frequency. This can 
be analyzed as a preferred mode in Ref. 10. Next, the 
perturbation increases the entrainment from the outer 
irrotational flow. This causes the reduction of curvature of 
the time-mean streamlines near separation line and it 
deduced the large pressure gradient Thus it decreases Cp, 
and C*. It is also found that the maximum of entrainment 
rate is observed at St,-. 

Many experimental results reveal that the effect of 
pulsation on flow structures is observed mostly at the 
comparatively lower frequency region 1,31112 . At a 
particular frequency, the vortex coalescence is 
surprisingly enhanced. It causes that the reduction of 
recirculating separation bubble and flow drag of the body. 

The distributions of the time-mean velocities U and 
V are displayed in Fig. 8. The difference between non- 
pulsating flow and pulsating flow is seen to be negligible 
near the separation line where vortices are just formed out 
However, the effect of perturbation is evident in the 
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Fig. 8 Distribution of U/U, and V/U, along the 
downstream distance 
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Fig. 9 Distribution of r.m.s. u/U, and r.ms. v/U, along 
the downstream distance 

reattachment region. The relatively large value of V at the 
separation line is due to the abrupt decrease of the 
curvatures of streamlines. Fig. 9 shows the r.m.s of 
velocity fluctuations. The turbulence level of perturbed 
flow is found to be considerably higher than that of natural 

flow. 
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Fig. 11 (a) Distribution of discrete-vortices in 
unperturbed flow at t = 288 HAJ, 


Fig. 10 Distribution of Reynolds shear stress -u'v'/U, 2 


The distributions of the Reynolds stress are shown in 
Fig. 10. In perturbed flow, it is noteworthy that Reynolds 
stresses have negative values at x/H-0,4 and 8. Hussain 13 , 
in his study on the coherent structure of turbulent shear 
layer, proposed that Reynolds shear stress can have 
negative values of vortical structure in a certain 
configuration where vortex coalescence process is 
dominant. Furthermore, since perturbation strengthens the 
fluctuating components in shear layer, so the mean 
velocity may be affected by these fluctuating velocities. 
The characteristic wave length of vortices can be given as 
follows 12 ; 


H Uj * 


( 12 ) 


where Uj is the convection velocity of vortices and St,** is 
the most amplified shedding frequency, which will be 
referred later. Suppose the value of lyUj was 0.5 M , then 
L has the value of about 5H. This agrees that the locations 
of negative Reynolds stress are found to be nearly 
identical to the multiples of L. This relation reduces that 
negative Reynolds shear stress are brought about by the 
enhancement of vortex coalescence by pulsation and the 
positions are closely related with the characteristic wave 
length. From this, it can be predicted that the vortex 
merging is observed at every characteristic wave length in 
separation bubble. 

The patterns discrete-vortices in the separation 
bubble at a certain time(t=288H/U i ) for both cases are 
plotted in Fig. 1 1. It may be noted that large-scale vortices 
are more evident in perturbed flow. In addition, the rolling- 
up behavior of large-scale vortices, the merging of 
vortices and the vortex-shedding from separation bubble 
to downstream are also observed. It is noteworthy that the 
concentration of discrete-vortices occurs at each multiples 
of L. It gives a numerical validation for negative Reynolds 
shear stresses in separation bubble. 
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Fig. 11 (b) Distribution of discrete-vortices in 
perturbed flow at t = 288 H/U, 



Fig. 12 (a) Width of shear layer with time at leading- 
edge in natural flow 



Fig. 12 represents the width of shear at separation as a 
function of time. Approximately, this width may be 
regarded as a distance between the core vortex and the 
wall of body. In non-pulsating flow, as time advances, this 
width becomes constant However, in perturbed flow, at 
St^, it oscillates with the period of 10Ht/Ui, which 
corresponds to the forcing frequency St^-0.1. In 
perturbed flow, the nascent vortices interact with the 
surface from the start of the rolling-up motion, so this is 
related to the resultant amplification of disturbances. 

The power spectra, which as defined as eqn.(13), was 
obtained at the edge of the shear layer where u’/Uj=0.02 
and is shown in Ftg.13. As pointed out previously, the 
power spectra of u near reattachment region also discloses 
the broad peak value around St=0.1. Thus, this 
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Fig. 13 (a) Power spectrum of u’ (=E„.) near separation 
line at the edge of the shear layer (r.m.s. u/U,»*0.02) 



Fig. 13 (b) Power spectrum of u' (=E„.) near 
reattachment position at the edge of the shear layer 
(r.m.s. u/U^-0.02) 
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Fig. 14 Power spectrum of p' (=E p .) near reattachment 
position at the edge of the shear layer (r.m.s. u/U^0.02) 

frequency(St=0.1) can be considered as the most amplified 
shedding frequency of large-scale vortices, which is also 
verified from the power spectra of pressure at the same 
position(Fig.l4). 

d(fHAJi) = (u' 2 AJ? , P' 2 4 pU?) (13) 

In passing, it is notable that the present result tends to 
have appreciably the same value as the results of Mabey 15 . 
It gives the relation of St X r «0.7 with almost the same 
value of the present result. This trend has also been 
noticed in other experi m e ntal res ults 0 . It is stressed here 
that the most amplified shedding frequency gives strong 
effects on the reduction of reattachment length, drag and 
enhancement of vortex coalescence in the separation 

bubble. 

The power spectra near separation line shows the 

broad peak at relatively low frequency(Fig. 13(a)). This is 
regarded as an evidence of the flapping motion of 
separation bubble, which seems to be attributed to 
unsteadiness l - 3 - 14 . Near reattachment region, the low- 
frequency peak is gradually diminished and the broad 
peak is centered upon around St**«0.1. Roshko 3 called 
this frequency as the shedding-type instability, which is 
the dominant frequency in separation bubble. 

In this study, we assumed that St** is identical to 

St This is attributed to the analogy as follows; If the 

pulsating frequency is higher than the initial Kelvin- 
Helmholtz frequency St** there is no instability region. In 
this sense, the slow disturbance can not be amplified. Thus 
St** becomes dominant around the value of St^. 

The way to estimate St** has been tried through a lot 
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of experiments. Roshko 3 obtained the relation fh/U,»0.08, 
where h is the asymptotic height, U, is the velocity at 
separation line and f is the shedding frequency 
corresponding to St^. In the present case, we obtained 
the value of fh/U,=0.078, which is nearly identical to that 
of Roshko. U, can be obtained from the two-dimensional 
Bernoulli's equation. 

^=<l-C,.) ,n 04) 

It should be noted that the flow geometry of Roshko 
is different from the present experimental configuration, 
however, the result gives nearly same trends. Roshko's 
result was obtained from the frequency of von Karman 
vortex sheet in a blunt circular cylinder. Levi 17 also 
suggested the Strouhal Law, which defined as fd/U i "0.16. 
In this case, if d=2h substituted, his law gives the exacdy 
same value of our result. From the aforementioned 
comparison, the tentative summary can be suggested that 
has its origin in the large-scale structures of 
turbulence, regardless of flow geometries. In this 
connection, it is believed that St,,,* is approximated by 
St [hr 1 . Thus, this forcing frequency enhances the spreading 
rate of flow and vortex coalescence and minimizes the 
drag effectively. 
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